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ABSTRACT 

Aims. 60% of the A star members of the 12 Myr old /3 Pictoris moving group (BPMG) show significant excess emission in the 
mid-infrared, several million years after the proto-planetary disk is thought to disperse. Theoretical models suggest this peak may 
coincide with the formation of Pluto-sized planetesimals in the disk, stirring smaller bodies into collisional destruction. Here we 
present resolved mid-infrared imaging of the disk of r] Tel (AOV in the BPMG) and consider its implications for the state of planet 
formation in this system. 

Mefhods. 77 Tel was observed at 11.7 and 18.3//m using T-ReCS on Gemini South. The resulting images were compared to simple 
disk models to constrain the radial distribution of the emitting material. 

Results. The emission observed at 18.3/jm is shown to be significantly extended beyond the PSF along a position angle 8°. This is the 
first time dust emission has been resolved around rj Tel. Modelling indicates that the extension arises from an edge-on disk of radius 
0.5 arcsec (~ 24 AU). Combining the spatial constraints from the imaging with those from the spectral energy distribution shows that 
>50% of the 18yum emission comes from an unresolved dust component at ~4 AU. 

Conclusions. The radial structure of the 77 Tel debris disk is reminiscent of the Solar System, suggesting that this is a young Solar 
System analogue. For an age of 12Myr, both the radius and dust level of the extended cooler component are consistent with self- 
stirring models for a protoplanetary disk of 0.7 times minimum mass solar nebula. The origin of the hot dust component may arise in 
an asteroid belt undergoing collisional destruction, or in massive collisions in ongoing terrestrial planet formation. 
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1. Introduction surements put this system at a distance o f 48pc, giving a pro- 
j ected offset of 192 AU for the companion. iBackman & Parescd 

Debris disks around nearby stars are often discovered through ^^^^ identified this star as a debris disk candidate based 

^ the detection of emission beyond that predicted to come fromflie j^^g measurements which indicates that excess emission is 

H , stelkr photosphere in the infrared or sub-millimeter (see fctll t at 12, 25 and 60fim. The excess measurement at 12A.m 

. ^ , I2Q0M for a recent review). The detected emission arises from particularly rare ( Aumann & Probst 1991), as debris disks 

the absorption and subsequent re-emission of starlight at longer ^^^.^ ^^^jj ehai-acterised by dust at several 10s of AU 

wavelengths by dust particles in the disks. These dust grains ^j^^^j^ ^^^j ^^^^^^ ^^-^^ ^ 

are believed to be continuously reproduced by the collisional ^^^^ ^^^^^^^ ^-^^^ ^^^^ ^nA^^ion resolved at thermal wavelengths 

destruction and/or sublimation of larger bodies (planetesimals), j^^^^ ^j^,^^ extended on 100s of AU scales {fi Pic, HR4796A, 

as the lifetimes of such small grains due to removal by radia- 49 ^^^^^ ^4^59^ 32297 Backman & Paresce 199| 

tion pressure or Poynting-Robertson (P-R) drag are small com- j^yawardhana et al. 1998; Wahhai et al. 2007; F isher etalfeOOOt 

pai-ed to the age of the star (Backman & Paresce 1993). Studying Moerchen et al.„2007a) . Nevertheless, examples of A stars with 

these debris disks can give clues about how the circumstellar dust confirmed to he within lOAU are known (e.g., ^ Lep, 

region has evolved during the star s formation and subsequent |Moerchen et al.ll2007bl) . and there remains debate as to whether 

hfetime, how any planetary system around the star may have ^^^^ ^^^^\on c an originate in the s teady-state destruction of 

formeiand_evolved and the current status of the system (see pjanetesimal belts (IWvatt et al.ll2007bl) or in collisions between 

^ " ^ " o , j^ • , , r , T,- • growing protoplan ets that are expecte d during terre strial planet 

77 Tel (HD 181296) is an AOVn member of the p Pictons forj^ation (Kenvo n & Bromk^lM) . IChen et all (2006) pre- 

movmg group, with an age estimated to be 12+^ Myr (age es- ^^^^^^ j^.^ observations of 77 Tel showing that the excess emis- 

timated for moving group in Zucker man et al. || 2001|) . It has an ^j^^ ^^^^^^ ^ ^i^gj^ ^mperature of dust as the spectrum 

M7/8V companion at an offset of 4 at a position angle of 160° ^hey interpret this as evidence for two temperatures 

to the star (Lowrance et al. | 2000|) . Hipparcos hsted parallax mea- ^^^^^ ^^^^ emission spectrum alone it is not possi- 
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ble to tell if this arose from dust in two separate radial locations 
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around the star, or alternatively two different grain populations at 
the same radial location. The aim of this paper is to use resolved 
imaging of the disk to allow a differentiation between these al- 
ternatives. In addition these images will confirm that the excess 
emission is in fact centred on the primary star and not the binary, 
although due to the brightness of the 12/im flux it is unlikely that 
the emission arises from dust around the M7/8 companion. 

The age of this system and membership of the moving group 
make it very interesting in terms of evo l utiona ry studies of de- 
bris disks. Recent work by ICurrie et"aLl (l2007b has shown that 
excess emission at 24yum around A stars increases from 5-10 
Myr, peaks around 10-15 Myr and then declines with age. Th e 
cause of the peak is not yet known (see review in Wyatil2008l) . 
but is suggested to be the delayed formation of Pluto-sized bod- 
ies further from the star, as these bodies are necessary to stir the 
disk so that c ollisions are destruct ive thus releasing large quan- 
tities of dust dKenvon & Bro mlev 2004). The 12 Myr yS Pictoris 
moving group has 5 A star members, 3 of which are known to 
have debris disks that exhibit exce ss emission in the nud-infrared 
06 Pictoris, HD 172555 and ij Tek lChen et al.ll2006HRebull et all 
[2OO8), and the proximity of this group means it offers the chance 
to assess the origin and diversity of the peak. 

Here we present 11.7 and 18.3//m imaging of rj Tel taken 
with the T-ReCS instrument on Gemini South which shows a 
resolved debris disk at 18.3/vm. We confront these observations 
with models to determine the optimal disk parameters and dis- 
cuss the implications of the inferred structure for the status of 
planet formation in t his system, both by compa rison with the 
self-stirred models of iKenvon & BromlevI (|2004|) and with our 
understanding of the early Solar System. 

2. Observations 

The source was observed under proposal GS-2007A-Q-45 us- 
ing T-ReCS on Gemini South with filter Qa (Ac - 18.3 jum, AA 
= 1.51 fj.m) and Si-5 (Ac = 11.66 fim, AA = 2.13 fJ.m), which 
we refer to as Qa and Si-5 observations henceforth. The pixel 
scale of the T-ReCS imager is 0'.'09 with a total field of view 
of 29" X 22". The observations were taken in parallel chop-nod 
mode with a chop throw of 10"and chop position angle of 55° 
(East of North) at Q, 100° at N. The Qa observations were taken 
over two consecutive nights ( P' and 2"'' July 2007), with the Si-5 
observations taken on one night only (12* July 2007). Total inte- 
gration time was 9120s (4560s on source) at Qa, and 912s (456s 
on source) at Si-5. Observations of standard stars HD 196171 
(spectra l type KOIII) and H P 179886 (spectral type K3III), both 
listed in lCohen et all (Il999l) . were made to calibrate the photom- 
etry and also determine the PSF of the observations. The order 
of observations and integration times are shown in Table [T] 

These data for rj Tel were taken in multiples of 304s-long in- 
tegrations, each consisting of 7 complete nod cycles (ABBA). 
The data were re duced using custom routines described in 
I Smith et aP (l2008t) . Data reduction involved determination of 
a gain map using the mean values of each frame to determine 
pixel responsivity (masking off pixels on which emission from 
the source could fall, equivalent to a sky flat). In addition a dc- 
offset was determined by calculating the mean pixel values in 
columns and rows (excluding pixels on which source emission 
was detected) and this was subtracted from the final image to 
ensure a flat background. Pixels showing high or low gain, or 
those which showed great variation throughout the observation 
were masked off. In order to avoid errors in co-adding the data 
which could arise from incorrect alignment of the images, we 
used re-binned images in which each pixel was a fifth the size of 



the T-ReCS image pixels to determine an accurate centre of the 
stellar images, aligned and co-added these images and then re- 
binned to the original pixel size. The re-binning was performed 
using bilinear interpolation across the array. This technique was 
used both to determine a final image of rj Tel and standard star 
for each band. To ensure that background noise did not affect 
the re-binning process the same procedure was applied with an 
additional step of smoothing the images with a Gaussian kernel 
of FWHM equal to that of the standard star image before re- 
binning. It was found that the relative sizes of the final images 
of 77 Tel and the standard stars and ellipticity of the images was 
unchanged by this additional step, indicating that the re-binning 
was not affected by background noise. This is as to be expected 
in data such as this where the detector super-samples the PSF. 



Date 


Object 


Filter 


Integration 


Calibrated 








time (s) 


flux, mJy 


01/07/2007 


HD 196171 


Qa 


152 


6058 


01/07/2007 


j/Tel 


Qa 


1824 


350+ 11 


01/07/2007 


HD 196171 


Qa 


152 


6058 


01/07/2007 


j/Tel 


Qa 


1824 


347 ± 1 1 


01/07/2007 


HD 196171 


Qa 


152 


6058 


02/07/2007 


HD 196171 


Qa 


152 


6058 


02/07/2007 


?/Tel 


Qa 


1824 


345 + 10 


02/07/2007 


HD 196171 


Qa 


152 


6058 


02/07/2007 


j/Tel 


Qa 


1824 


350 ± 11 


02/07/2007 


HD 196171 


Qa 


152 


6058 


02/07/2007 


j/Tel 


Qa 


1824 


331 ± 10 


02/07/2007 


HD 196171 


Qa 


152 


6058 


12/07/2007 


HD 179886 


Si-5 


304 


4417 


12/07/2007 


?/Tel 


Si-5 


912 


342 ± 18 


12/07/2007 


HD 179886 


Si-5 


304 


4417 



Table 1. Observations in order and total integration times of the 
observations taken under proposal GS-2007A-Q-45. Note that 
on-source integration time is half the total integration time listed 
in the table. 



3. Results 

The final co-added images of rj Tel in both filters and the stan- 
dard star for PSF reference are shown in Figure [1] The con- 
tours shown in the images are at 20%, 40%, 60% and 80% of 
the peak in each case (peak value in Qa 5.24mJy/pixel, in Si- 
5 9.89mJy/pixel). Also shown in the right-hand image for each 
row is the residual emission around rj Tel observed when sub- 
tracting the standard star image as a PSF reference scaled to the 
peak of the 77 Tel image. The contours and residuals are discussed 
further below. 

Photometry was performed using I'.'O radius circular aper- 
tures centred on the stellar image. The centre of the image was 
determined through a 2-dimensional Gaussian fit. An average of 
the calibration factors determined from the standard star obser- 
vations taken on the night was used for corresponding photomet- 
ric calibration. Statistical noise was determined using an annulus 
with inner radius 1 "and outer radius 2"centred on the source, re- 
sulting in uncertainties of 0.13 and 0.06 mJy/pixel at Qa and Si-5 
respectively. The calibration factors were found to vary by 4% 
and 2% on the first and second nights of Qa band observations, 
respectively, which is consistent with the variation of 3% found 
between fluxes measured on individual 304s integrations on 77 
Tel (which is bright enough for an ~60 sigma detection even in 
such a short integration). Although the calibration factors varied 




Fig. 1. The co-added final images from the Qa observations (top row) and the Si-5 observations (bottom). Orientation of the images 
is North up. East left. Left: The total co-added images of rj Tel. Middle: The co-added images of the standard star forming our 
reference PSF. Note the distinct ellipticity with major axis in the North-South direction seen in the Qa image of rj Tel is not seen in 
Si-5 or in the standard star images. Solid line contours are at the 20%, 40%, 60%, and 80% of the peak levels for all images. Peaks 
are at a level of 650 mJy/arcsec^ for rj Tel at Qa, 1220 mJy/arcsec^ for rj Tel at Si-5, 13990 mJy/arcsec^ for the standard star at Qa 
and 20670 mJy/arcsec^ for the standard at Si-5. Right: The residual emission after subtraction of the total standard star image scaled 
to the peak of the image of rj Tel. The dashed contour indicates the 3cr per pixel level (from background noise). The colour scale is 
linear from to the maximum pixel value in each image. The residual images have peak brightness of 88 mJy/arcsec^ at Qa and 35 
mJy/arcsec^ at Si-5. The arrow indicates the direction of the M7/8 V binary which lies at 4", PA° 160. Note that the emission seen in 
the Si-5 residual image is not significant (< Icr significance) when the variation in the PSF is considered. 



by just 5% for the two observations of standards in the Si-5 fil- 
ter, analysis of the individual 304s integrations on 77 Tel (each of 
which resulted in a 78 sigma detection) showed that photometric 
accuracy was in fact at a level of 13%. Our photometry yields to- 
tal fluxes of 342 + 44 mJy in Si-5 (with S/N of 292) and 345 + 15 
mJy in Qa (with S/N of 129), including both calibration and sta- 
tistical noise, where the S/N given in brackets indicates the level 
of statistical noise. Standard stars were observed at a similar air- 
mass to 77 Tel. In addition, calibration levels were compared to 
airmass for the standard star observations and no correlation was 
found, therefore no correction for extinction was applied to cal- 
ibration. Using a Kurucz model profile scaled to the 2MASS K 
band magnitude, the stellar flux in filters Si-5 and Qa is expected 
to be 282 mJy and 1 14 mJy, respectively. Thus the photometry 
gives an excess of 60 + 44mJy and 23 1 + 15mJy in Si-5 and Qa. 
The IRS spectrum of ' Chen et al.l (l2006l) after subtraction of the 
above mentioned photospheric emission gives excess emission 
of 150+5 and 257+2 mJy, respectively, at these two wavelengths 
and so the photometry presented here is in agreement within the 
T-ReCS photometric errors with the IRS results at the 2cr level. 
Thus it does not appear that the IRS spectrum includes any emis- 
sion centred on the star outside the 1 "aperture used here for pho- 



tometry, as might happen in the larger IRS beam (extracted along 
a 3'.'7 sHt). 

Figure [T] shows the final co-added images of 77 Tel and the 
standard star for a PSF reference. The top line of this figure 
shows the Qa total image of rj Tel and the standard. An elliptical 
shape can clearly be discerned, and a 2-dimensional Gaussian 
fit to the image has an ellipticity of 0.1 17 + 0.007 with the ma- 
jor axis at 19.9 + 2.0° East of North; the same ellipticity was 
seen on the two separate nights of Qa observing. For compari- 
son, the ellipticity of the final standard image is 0.084 + 0.005 
with the major axis at 83+2°. This slightly elliptical shape of the 
PSF is seen across the standard star integrations (see later this 
section) and does not coincide with the extension seen in the 77 
Tel images. As all Qa observations were performed with no on- 
sky rotation and chop-nod performed at 55°, the direction of the 
extension is unlikely to be the result of chop smearing. 

Also shown in Figure[r|is the residual emission of 77 Tel fol- 
lowing a subtraction of the PSF reference standard star image 
scaled to the peak of the 77 Tel image, which subtracts a total 
of 252mJy from the image. Although this subtraction method 
is likely to remove some of the disk emission (some of which 
contributes to the peak surface brightness), it does highlight the 
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-1,5 -1.0 -0.5 0.0 0.5 1,0 1.5 -1.5 -1,0 -0.5 0.0 0.5 1.0 1.5 

Arcsec Arcsec 

Fig. 2. The profiles of line cuts through the total images, at 8° (left) and 98° (right) at Qa. The median FWHM measurements taken 
over all sub-integrations from Gaussian fits to the profiles at the two angles are indicated for reference. Also shown are averaged 
profiles from all the line cuts for both rj Tel and the standard star. Observe that the profile for Tel at 8° is much broader in the 
wings than the standard star The Icr error bar shows the size of +lcr background error per pixel from the rj Tel image. Note that 
multiple intensity profiles for the science target and standard are shown here. These are taken from longer integrations than noted 
in the text to reduce the number of profiles and improve the clarity of the figure. 



location of the non-pointlike emission. Peaks are seen in the Qa 
residual image that are centred at 0'.'45 (21.6AU), 5.2° East of 
North, and a.'40 (19.2AU), 190.5° East of North. Symmetrical 
peaks in the residual emission either side of the star are exactly 
what would be expected for an edge-on ring dTelesco et al.ll2000l 
although also note that we cannot rule out a less inclined ring 
without detailed modelling of the disk, see section 4.). Since the 
residual peaks are more representative of the extension seen than 
the 2-dimensional Gaussian fit discussed above, hereafter the ex- 
tension is considered to be at a position angle of ~8°. Thus the 
extension is at ~28° with respect to the direction of the binary 
(shown with an arrow), suggesting if this arises from an edge-on 
disk, then that disk is not in the same plane as the binary. The 
flux in the residual peaks totals 80 + 5 mJy (measured in 0'.'4 
radius apertures; the error term includes background and cali- 
bration uncertainty added in quadrature), which as expected is 
below that of the excess emission 231 + 15mJy measured in our 
photometry. An estimate of the level of disk emission which may 
have been removed from the residuals image by our method of 
PSF subtraction can be determined from consideration of how 
much the two peaks may have contributed to the central pixel. 
Using the PSF model scaled to the residual peaks this was calcu- 
lated as 0.3mJy (as compared to the total peak before subtraction 
of 5.24mJy). This method suggests the extended emission repre- 
sents a minimum of 6% of the peak emission, as if the extended 
emission is truly a disk lying across the star, this method will 
again underestimate the disk flux in the central pixel). A more 
detailed model of the extended emission seen in the Qa image is 
discussed in section 4. 

To assess the temporal variability of this FWHM and the sig- 
nificance of the extensions, as wefl as to search for a dependence 
on time or airmass, linecuts were performed at position angles 
of 8° and 98° for individual nod set images (22s) for the stan- 
dard star and two nod set images (43s) for 77 Tel. Two methods 
of fitting the cuts, by a Gaussian profile and a Moffat profile 
were tried. Both fits provided the same relative results, although 
the peak was normally fitted more accurately using the Gaussian 
profile. No dependence on airmass or time was found, and me- 
dian FWHM values and standard errors were found to be: for the 
standard target 0'.'531 ± 0'.'002 at 8°, 0'.'566 ± 0'.'004 at 98°; and 



for 77 Tel a.'582 + a.'OlO at 8° and a.'545 + tt'OO? at 98°. The 
lower signal-to-noise on the 77 Tel observations naturally leads 
to a larger degree of scatter in the intensity profiles. However 
it is clear that at Qa there is no difference in the FWHM of the 
77 Tel frames and standard star target frames at 98°, but that at 
8° the FWHM is larger for 77 Tel than for any of the individual 
integrations of the standard stars. A further graphical demonstra- 
tion of the significance of the extension is provided by Figure |2] 
which shows the cuts through each Qa integration (of 304s) from 
the standard star imaging (plus signs) and the 77 Tel observations 
(dots). At 8° the 77 Tel sub-integrations are not only broader at the 
FWHM but also have broader wings to their profiles compared 
to the standard star images. Neither of these features is seen at 
98°. 

The Si-5 image of 77 Tel shown in Figure [1] includes only 
the first 2 of 3 integrations on the source. The final integration 
showed an unusual shape compared to the previous integrations: 
a greater extension in the North-South direction, and an unusual 
asymmetric shape, with an extension towards the East of the im- 
age not matched on the West side. This shape was seen in all 
but the first nod set. Neither the North-South extension nor the 
Easterly extension match the chop PA of 55°. The eccentrici- 
ties and orientations of the major axis for the 3 integrations are: 
0.027+ 0.009 at 1 10°; 0.082 ± 0.009 at 167°; and 0.163 + 0.012 
at 80°. Thus we consider that the shape seen in the final integra- 
tion is not evidence for true extension, but is due to the variable 
seeing conditions on this night. 

In the Si-5 residual image we see only low significance emis- 
sion at the Northern edge of the original stellar image. The emis- 
sion seen in this residual image (Figure [T] lower right image) is 
6.6 + 0.9 mJy (this error includes only background error, not 
PSF uncertainty), however with subtraction of individual stan- 
dard star images rather than the coadded standard star PSF refer- 
ence image flux in this region can be as low as -7.8mJy, demon- 
strating the importance of consideration of the PSF variation 
when determining evidence of extended emission close to the 
star, and the true insignificance of this emission. Again the line- 
cuts for individual nodsets, both for 77 Tel and standard, were 
assessed. No dependence on airmass or time was found, and me- 
dian values and standard errors of the Gaussian FWHM mea- 
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surements excluding the final integration for rj Tel (see above) 
are: for the standard target a.'383 + (T/OOS at 8°, 0'.'407 + Of.'OQG 
at 98°; and for // Tel a.'386 + a.'003 at 8° and 0'.'414 ± Qf.'OU at 
98°. Thus examination of the PSF uncertainty indicates that the 
residual emission is likely to be the result of PSF variation over 
the course of the observation. 

4. Modelling the emission 

In order to determine the constraints on the radial distribution of 
emission provided by the resolved Qa imaging, we considered 
models for disk structure of increasing complexity. For all mod- 
els it is assumed that the images are comprised of unresolved 
stellar flux (of 114 mJy), plus an additional unresolved disk com- 
ponent (motivated by the hot component inferred by Chen et al. 
2006) of flux Funres, along with an additional resolved disk com- 
ponent of flux Fres- The IRS photometry was used to constrain 
the total flux so that F^nres + ^^res =257mJy. The resolved flux 
is assumed to arise from an axisymmetric disk of opening angle 
5° (which is not constrained by the modelling) which lies at an 
inclination / to our line-of-sight. Since the emission is consistent 
with an edge-on disk / = is assumed in the modelling, although 
the extent to which it is edge-on is constrained later (see section 
4.2). Two types of resolved disk were assumed. Narrow rings 
are defined by r, the radius of the mid-point of the disk, in AU; a 
finite ring width of dr/r - 0.2 was assumed, with a constant sur- 
face density. Extended rings were defined by a continuous distri- 
bution of surface density from some inner radius to an outer 
radius rout with surface density following a power-law E oc r'', 
assuming grains with black body temperatures and emissivities. 

The total models (disks and point sources) were convolved 
with the average PSF (shown in Figure [T]) for comparison with 
the observed data. The compatibility of each model with the ob- 
served structure was tested by taking line cuts along the direction 
of extension (at a PA of 8°). The line cuts were averaged over 3 
pixels (0'.'27) in width centred on the peak of the emission, taken 
over a length of 50 pixels (i.e., + 2'.'25 from the centre) to give a 1 
dimensional profile for both the observation of rj Tel (with error 
bars) and for each model tested (e.g., Figure[3] right). For each 
model described below a range of parameters was tried and a 
goodness-of-fit test applied to find the best fitting parameters 
and confidence limits (errors from background statistical error). 
In addition to the from the line-cuts, that from a comparison 
of the model images and observation was also considered (e.g.. 
Figured left). 

To test whether more complicated models are necessary to 
explain the images, a Bayesian Information Criterion (BIC) test 
was applied to the model fitting: 

BIC^Nln(xij + k\nN, (1) 

where is the number of data points and k the number of 
free parameters, and xl,:„ is the minimum value fou nd from 
each model (Schwarz i978; see also Liddle 2004; Wahha Tet alj 
l2005h . This test penalises models for having unnecessary param- 
eters. A lower value of the BIC is preferred; a difference of 2 
between BICs of different models indicates positive evidence 
against the higher BIC value, and a difference of 6 indicates 
strong evidence against the higher BIC value. 

4.1. One parameter models 

The simplest model is one with no resolved disk component. In 
this case we allowed Fumes to be a free parameter, but as expected 



from the images and comparison of radial profiles in Figure |2] 
the testing showed that no model provided a good fit to the 
profile with ;t'^j^=786 and a reduced ;^f^=16, where is a 

perfect fit. The BIC of this model was 98 (from the 1 dimensional 
profiles). 

Another possibility is to assume that there is no unresolved 
disk component, just a narrow ring defined by the free parameter 
r. Values of r tested in this case were in the range 8-60AU. This 
lower limit was chosen as smaller disks appear unresolved and 
thus provide the same goodness-of-fit as the unresolved emission 
only model described above. In this case the fit to the observed 
profile was poor, worse than assuming the emission to be com- 
pletely unresolved with ;t'^j^=2815 and a reduced The 
BIC of this model was 160. This demonstrates that the images 
confirm the suggestion of Chen et al. that the material is at mul- 
tiple radii. 



4.2. Two parameter (preferred) model 

The two free parameter model assumes that the disk emission is 
comprised of a narrow ring of radius r and an unresolved com- 
ponent at a level Fumes (and so is essentially a combination of 
the two one parameter models considered above). The free pa- 
rameter Fumes, which also controls the value of Fres, was varied 
between and 280 mJy (in 35 mJy steps), while r was tested 
with values between and 95AU (in 3 AU intervals). A much 
better fit was achieved with the two parameter model with an 
unreduced y^. =53 and a reduced =1.13 for r = 24 AU 

^ mm ^ rmm 

and Fumes = 205mJy (corresponding to F^es = 54mJy). The pro- 
file of this best fitting model is shown in Figure |3] in which the 
contributions of the different components are identified. The BIC 
for this model was 31 (from 1 dimensional profiles) indicating 
that this is strongly favoured over the 1 parameter models. 

The x^ minimisation showed that the parameters were con- 
strained to the ranges: Fumes = 116 - 225, 136 - 231, 96 - 243 
mJyandr = 16-32, 10-42, 10-53AU at the 1,2, and 3o- levels 
respectively (ranges from a marginalisation of the joint distribu- 
tion of the two parameters). Using the line cuts perpendicular to 
the direction of maximum extension the inclination to our line 
of sight was constrained to be less than 20°, since at / > 20° the 
model FWHM (a.'574) exceeds the observed FWHM (0'.'566 + 
0'.'004) by more than Icr. 

Given the importance of PSF variation in the residual struc- 
ture seen in the N band images, and as we used the average 
PSF in the modelling outlined above, we performed the same 
modelling process using the different PSFs measured in the sin- 
gle standard star observations (i.e. not co-added). The sensitiv- 
ity of our best fit parameters on the chosen PSF was minimal, 
with typical 1 cr confidence regions close to those given in sec- 
tion 4.2 (largest ranges were 8AU < r < 32AU and 203mJy 
< ^^umes < 339mJy). Thus the variation in the observed PSF does 
not strongly change the best fitting model for the extended emis- 
sion. 

The final total model at Qa is shown after subtraction of the 
PSF scaled to the peak of the observed image of rj Tel in Figure[3] 
along with the profile of this best fitting model with an edge-on 
disk (inclination = 0). Subtracting the model image from the ob- 
servation reveals that two are a close match within the statistical 
noise over the area of the image (900 pixels = 2'.'7x2'.'7, reduced 
X^ = 1.3). 
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Fig. 3. The best fitting two parameter model for the Qa band emission around rj Tel. This model includes that stellar flux, an 
unresolved flux component at 205 mJy and and an edge-on resolved disk component between 21-26AU. Left: The model "residual" 
image, with contours from the observed residuals image (Figure[T]top right). The model was convolved with the PSF and subjected 
to the same subtraction (standard star PSF reference image scaled to peak) as the rj Tel image. The scale is the same as the observed 
residuals image. The observed residuals contours show the good correspondence of the model with the observed emission and are at 
1/3 and 2/3 of the residual peak. Right: The line cut through the rj Tel Qa image at 8° and the best fitting model at the same position 
angle. 



4.3. More complicated models 

As soon as the ring is assumed to be broad the number of free 
parameters increases to 3 or 4, depending on whether an unre- 
solved component was included: i.e., ^n, rout, and possibly 
Fumes- We Started by assessing whether the radial distribution 
could be explained by a single continuous distribution of dust, 
without an unresolved component. We tried model grids with 
Hn e (0.2, 6)AU spaced at 0.2AU intervals, rout e (6, 60)AU in 
1 .5 AU intervals and y tested at intervals of 0.5 with y e (-2, +2). 
This resulted in a poor fit to the data (^^=198 over the 1 dimen- 
sional profile, reduced ;t'^=4. 2), and the BIC of 83 which sug- 
gests that the radial distribution is discontinuous. 

It was possible to find a better fit to the observations than 
the two parameter model of section 4.2 when F^nies was also 
included as a free parameter; a reduced of 1.03 (from the Id 
profile) was found for ri,, = 1.5;^q5AU, rout = 45 ± 5AU, y - 
-0.5 and Funre.s = 170+26 mJy. However, the BIC calculation 
results in a value of 55 suggesting that the additional parameters 
required for an extended dust distribution are not justified to fit 
the data. 



4.4. Combining imaging and SED constraints 

The emission spectrum of rj Tel is shown in Figure |4] Here the 
stellar photospheric emission has been modelled using a Kurucz 
atmosphere of 9506K that has been scaled to the 2MASS K band 
photometry (implying a stellar luminosity of 22L0 for 48pc). 
The photosphere of the M7.5 star companion is shown for ref- 
erence and is modelled using 2700K NEXTGEN model atmo- 
sphere (iHauschildt et al.l [19991) . This companion is outside the 
aperture used to determine the TReCS photometry. In any case 
the emission from this M7.5 star at the wavelength range cov- 
ered by both TReCS and the Spitzer IRS spectrum is below the 
background statistical noise level and thus it does not contribute 
to the photometry shown at >5yL(m. Photospheric contributions 
have been subtracted from the fluxes shown at > 5yum which 



represent the excess emission from the disk. The results of our 
imaging require that the emi tting dust be at two radii, in agree- 
ment with the conclusions of lChen et alj (l2006l) . These separate 
components are shown by open symbols on Figure H] 

Modelling was performed to provide a two component fit to 
the SED in which each component (the unresolved excess and 
resolved disk) is pegged to its 18yum flux, with the unresolved 
component dominating the IRS spectrum < 18;um and the re- 
solved component dominating the longer wavelength fluxes. We 
assume a maximum grain size of lOOOjum, with a minimum grain 
size of 10.8yL(m as this is the size at which grains are removed 
from the system by radiation pressure. Large grains (greater 
than mm-sized) are inefficient absorbers and emitters of radia- 
tion, and their addition to the grain population would not have a 
great affect on the SED. A size distribution of n(D) oc D^^ ^ (see 
Dohnanyi 1969) was also assumed. We consider a distribution of 
grain sizes to be a realistic model for a dust populatio n created in 
the different origin scenarios presented in section 5 dChen et al] 
2006 use a single grain size to fit the IRS spectrum). Grains 
were assumed to have a silicate core and organic refractory/icy 
mantle (as used in Li & Greenberg 1998; Aug ereau e t al. 1999i) . 
Dielectric constants were calculated using Maxwell-Garnett ef- 
fective medium theory and optical properties calculated using 
Mie theory, Rayleigh-Ga ns theory and Geometric O ptics in the 
appropriate size regimes dBohren & Huffman|[T983l) . The unre- 
solved component was then best fit by a silicate fraction of 30% 
by volume (with the remaining volume composed of organic re- 
fractory material) and porosity of 0.2, which fit the shape and 
level of the silicate feature at ~ lOjum, resulting in a dust popu- 
lation at 3.9AU of fractional luminosity Ldust/^-* = 1.57 x lO"'*. 
Although we do not consider this modelling to have uniquely 
constrained the grain composition, this is a physically based 
model from which it is possible to estimate the radial location of 
the dust. Since this location is smaller than the FWHM it is con- 
firmed that this should appear unresolved in our images, which 
would be expected as long as the radius was < 6AU. 
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Fig. 4. Spectral energy distribution (SED) of 
J] Tel. The photosphere of 77 Tel is fitted with 
a Kurucz model profile scaled to the 2MASS 
K band flux and shown with a solid line, as 
is the photosphere of the binary companion 
which is modelled with a NextGen model atmo- 
sphere (appears on this plot only shortwards of 
5/im). Observations at > 5//m are shown after 
the subtraction of p hotospheri c contributions: the 
IRS spectrum of Chen et alj (120061) (sohd line 
between ~5-30//m), MIPS fluxes dRebuU et alJ 
I2OO8I crosses), IRAS colour-corrected excess (as- 
terisks) and the contribution to the 18yum fluxes 
from the unresolved (open triangle) and resolved 
(open diamond) components, as inferred from the 
modelling of section 4.4 (see also Table |2]i. Note 
that the IRS excess emission spectrum at < 8//m 
is subject to large uncertainty due to uncertain- 
ties in the level of photospheric emission in this 
range. The two component fit to the excess spec- 
trum is shown by dotted lines - the hotter unre- 
solved component at 3.9AU and a resolved com- 
ponent at 24AU. 
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Table 2. The observed emission from the T-ReCS observations of r] Tel, with comparison to the predicted emission from the two 
components used to fit the excess SED as described in section 4.4. Recall that the Si5 photometry suffered from poor calibration 
accuracy, and as such SED modelling was primarily fitted to the IRS spectrum at shorter wavelengths. Figure|4]shows the fit recreates 
the IRS spectrum accurately. 
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The same dust grain size distribution was used to model the 
cold resolved component. The same dust composition with the 
dust located between 21 and 26 AU as constrained by modelling 
the resolved Qa emission provides a good fit to the Spitzer IRS 
spectrum of the excess emission (Figure HJi. The fractional lu- 
minosity of the resolved component is similar to the unresolved 
component at Ljust/i'* = 1.39 x 10 The N band flux of the re- 
solved component is just 6mJy, which is consistent with the non- 
detection of the resolved emission in the imaging which is dom- 
inated by the unresolved component (at 140mJy). The full SED 
fit was compared to the observed photometry and ax^ goodness- 
of-fit computed. Data points from IRAS, MIPS, TReCS and the 
IRS spectrum were weighted according to their uncertainty. Data 
points at <7//m were excluded from the fit due to high levels of 
uncertainty arising from the correct determination of the stellar 
contribution at this range. The resulting reduced chf of 1.14 in- 
dicates the model fits the observed emission well. 



Thus the two component model provides a consistent pic- 
ture of this disk (see Table |2] and Figure |4] for comparison be- 
tween model and observed emission levels), and while we can- 
not accurately assess the range of temperatures (and radii) that 
are present in the two components, the poor fit of a continuous 
surface density distribution (section 4.3) suggests that there is a 
region of lower surface density between the two. 



5. Discussion 

In the following section we explore the possible origin of the two 
dust populations around t] Tel. We consider two possibilities: that 
7] Tel represents an analogue of the Solar System when it was at 
a similar age; and that the observed emission could arise from 
ongoing planet formation. We discuss the implications of these 
scenarios for any as yet undiscovered planets in the system, and 
place the 77 Tel disks in a wider context by comparison to the 
other/? Pic moving group A star disks. 



5.1. A young solar system analogue 

The two component model for the rj Tel disk represents an in- 
triguing analogue to the young solar system. Current theory sug- 
gests that the asteroid and Kuiper belts in the Solar System were 
originally much more massive than they are today and that much 
of the mass depletion occurred during a chaotic period known 
as the Late Heavy Bombardment several hundreds of Myr after 
the formation of the Solar System, possibly triggered by the mi- 
gration of the gi ant planet s (Gomes et al. 2005). According to 
the Nice model (iTsiganis et al.l 120051) . which reproduces fairly 
accurately the current confi guration of the Sola r Systems giant 
planets and the Kuiper Belt ( Levison et al.l200'7l) . the Kuiper belt 
was situated between ~15.5-34AU, with the orbits of the outer 
planets confined to the 5-15AU region, whereas the asteroid belt 
was at its current location. Thus in terms of radial location it 
is tempting to interpret the rj Tel disk in terms of an analogous 
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asteroid belt (the unresolved component at 3.9AU) and an anal- 
ogous Kuiperbelt (the resolved component at 24AU). It is worth 
noting that the luminosity of rj Tel is much higher than the Sun 
at ~ 24Lq, and thus the snow line (the location at which ices are 
expected to form) would be much further from the star. Using 
a simple 150K blackbody approximation the snow line would 
be expected at ~17AU around t] Tel (compared to 3.5AU in the 
Solar System). 

The current levels of excess emission arising from the two 
populations of debris in the Solar System are 2-3 orders of mag- 
nitude lower than that observed around 77 Tel, since the frac- 
tional lum inosity of asteroid belt dust is Ldust/^-* ~ 10"'^ - 10"^ 
jPermott et al. 2002), and that of Ku iper belt dust is Ldust/^-* ~ 
10-7 _ 10-6 dstern & Colwelllll997l) . However, the luminosities 
would have been much larger at an age similar to 77 Tel, since 
this would have been before the depletion of the asteroid and 
Kuiper belts. The original belts in the Nice model had masses of 
~lMe and ~35Me respectively (total mass in full size range of 
bodies), and so would also have had correspondingly higher lev- 
els of dust emission. The fractional luminosity of the Kuiper belt 
was estimated to be around 5 x 10 during the pre-LHB phase 
(Wyatt 2008; Booth et al. in prep.), which is comparable to that 
of the resolved component of the 77 Tel disk. 

5.2. Origin of fhe unresolved component 

The fractional luminosity of the pre-LHB asteroid belt would 
also have been considerably higher, since IM® of material at 
~ 3AU distributed in a collisional cascade size distribution 
from 2000km asteroids down to 0.8/vm dust blown out by radi- 
ation pressure would have resulted in a fractional luminosity of 
~ 10"^. However, collisional processing w ould ensure that such 
a high mass would not be long-lived, and IWvatt et alJ d2007bh 
showed how there is a maximum fractional luminosity for plan- 
etesimal belts of a given age and radius. For a belt at 3.9AU 
around rj Tel that has been undergoing catastrophic disruption for 
12Myr this maximum fractional luminosity is /max = 4 X 10"^. 
Although the current fractional luminosity is 30-40 times larger 
than this maximum value, this does not rule out the possibility 
that the unresolved component arises from the steady state col- 
lisional destruction of an asteroid belt for two reasons. First, the 
model requires /obs/ /max » 100 times larger to indicate a dis- 
crepancy given the various uncertainties in the model parameters 
including assumptions about the maximum planetesimal size, 
planetesimal strength and orbit al eccentricities o f around 0.05 
used to determine /„ax (see also Lohne et alj|200 8^. Second, it is 
possible that collisions in the asteroid belt only became destruc- 
tive after the protoplanetary disk dispersed (after which time ec- 
centricities can increase due to reduced dampening), and so more 
recently than 12Myr thus increasing /max- 

However, another possible interpretation of emission at a 
few AU is that it arises i n the late stages of the formatio n 
of terrestrial planets (e.g., Rhee et alj|200'7l: iMever et alj|2008h . 
High levels of hot emission are expected in the planet forma- 
tion models of iKenvon & BromlevI ( l2004l) once the planetesi- 
mal disks are stirred by the formation of Pluto-sized objects. 
As such objects form within 10^ yrs at 4AU in a quiescent 
disk with the same density as the m inimum mass solar nebula 
(MMSN) (iKenvon & BromlevI l2005h . stining of the disk by a 
newly formed Pluto does not seem likely to explain the hotter 
dust emission. However, evidence from the Solar System sug- 
gests massive collisions between terrestrial planet cores could 
occur quit e frequently in the first <100 Myr (e.g. formation of 
the Moon. ICanupll2004h . and models indicate that the final stage 



of ter restrial planet forma tion is characterised by such collisions 
(e.g., iKokubo et alj|200 ^. Thus the hot emission around 77 Tel 
could be evidence of stochastic collisions between large rocky 
cores. 

The unresolved component does not represent a population 
of PR-dragged particles from the resolved ring. Consideration 
of the collisional and PR-drag timescales for particles fr om the 
resolved ring (fcou/fpR ~ 0.05 for a belt at 3.9AU; see IWvati 
I2OO5L for equations) indicates that any grains susceptible to PR- 
drag released by collisions in the outer belt are likely to undergo 
further collisions before PR-drag can affect their orbits to move a 
significant dust population to the radius of the unresolved com- 
ponent. Furthermore, were PR-drag a dominant process in this 
disk system, we would expec t a dust distribution with a constant 
surface density (IWvattll2005h . not the two separate populations 
observed in the Qa images. 

It is not possible to distinguish between the asteroid belt and 
protoplanet collision possibilities at this stage, although detailed 
spectral modelling may be able to determine the dust composi- 
tion and so give clues as to the nature of the parent body (e.g., 
iLisse et aLll2007l) . Resolving the disk could also aid in its in- 
terpretation, although the small spatial scale of this population 
would make this a challenge. To date only the asteroid belt ana- 
logue aro und 4" Lep has been reso lved at a similar offset from the 
host star jMoerchen et al.ll2007bl) . 

5.3. Origin oftiie resolved component 

Although the suggested location and fractional luminosity of the 
resolved component around 77 Tel compare well with the early 
Kuiper belt, this does not require the two systems to have un- 
dergone a similar evolution and to have similar planetary sys- 
tems. An alternative scenario for the evolution of the 77 Tel 
syst em is give n by the self-stirred planet formation models of 
Kenyon & Bromley (2004) which may also explain the delay in 
the o nset of 24/zm emission peak until 10-15Myr jCurrie et al.l 
|2007 t). In such models planet formation within an extended plan- 
etesimal belt results in a bright ring of emission at the radius 
where Pluto-sized objects have recently formed. The timescale 
for the formation of Plutos, 600(E()/2mmsn) ' '"^ yr, around 77 Tel 
suggests that to have reached this stage at 24 AU by 12Myr re- 
quires a disk surface density equivalent to ~ 0.7 MMSN. 

To assess the possibility that the bright ring in the 77 Tel 
disk is caused by the recent formation of Pluto-sized objects at 
24 AU, we performed additional modelling of the 18/vm image 
based on the prescription for the evolution of dust surface density 
given in Wyatt (2008) t hat provides an emp irical fit to the self- 
stirred models of K envon & BromlevI ( 12005) . This model con- 
siders an extended planetesimal belt of surface density scaled to 
the MMSN in which annuli at different radii have a suppressed 
dust content until the age at which Plutos form at whi ch point 
their mass is eroded through steady state evolution (Wya tt et"al] 
|2007b). The dust emission is calculated assuming a ratio of dust 
area to planetesimal mass appropriate for a collisional cascade 
size distribution, and further assuming black body temperatures 
and emission properties. To account for possible global devia- 
tions from this size distribution or from pure black body grains, 
an additional scaling factor rj was applied to the surface bright- 
ness of the whole disk. An unresolved disk component was in- 
cluded in addition to the unresolved star The total flux from the 
self-stirred disc (as fixed by the disk surface density and scal- 
ing factor) and the unresolved component was fixed to 259 mJy 
from the IRS photometry (see section 4). Thus the model has 2 
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Fig. 5. The self stirred model of the resolved disk component 
as described in Section 5.3 after the same PSF subtraction as 
applied to the observations (Figure [T] Qa band residual image) 
and 2 parameter model, shown in Figure |3] The scaling in this 
image is the same as for Figure[T]top left. Contours are from the 
observed residuals image, at levels of 1/3 and 2/3 of the peak, as 
for Figure [3] 



free parameters; Zq/Smmsn, and r], with F^, 
to the level of flux in the self-stirred disk. 



fixed in response 



As expected from the simple calculation above, the best fit- 
ting model has a surface density of 0.7 times MMSN, leading to 
a ring at ~ 24 AU, with Funres = 216mJy (similar to the model 
of section 4.2) and 77 = 0.22. The low value of rj could arise 
because the grains are cooler than black body, since the ratio of 
By(l05K)/By{l23K) is ~ 0.3 at 18;um. The for this model 
is 58 with a reduced of ;if^j^=1.21, slightly higher than that 
of the 2 parameter model (x^^^^-l.l3) in Section 4.2. The re- 
duced of the image (taken as described in Section 4.2) was 
1.15, better than the two parameter model which had;t'^^jj^=1.3. 
The BIC for this model is 32, which is indicative of the almost 
equally good fit of the two models. The BIC values show there 
is no evidence against either the simple ring plus unresolved flux 
model (section 4.2) or the self-stirred model. The residuals for 
the best fitting self-stirred model are plotted in Figure |5] Multi- 
wavelength imaging may help break this degeneracy. The self- 
stirred model requires an extended disk with a brightness peak 
at the location of recent Pluto formation. If a larger disk location 
were observed in longer wavelength imaging this could be evi- 
dence in favour of the self-stiiTed interpretation, whereas in the 
solar system analogy we may expect to see the disk confined to 
a narrow radial range. Alternatively searching for planets in the 
system could further test these models. 

In other words the structure of the disk is entirely consistent 
with a self-stirred extended planetesimal belt in which the inner 
region has been cleared through collisional erosion following the 
formation of Pluto-sized objects. One caveat is that an additional 
unresolved component still has to be invoked. To be consistent 
such a component could not originate in an asteroid belt (unless 
this had significantly different material properties to the outer 
regions), and further work is needed to assess the duration and 
radial location of any ISyum emission expected from the final 
collisional phase of terrestrial planet formation. 



5.4. Implications for the planetary system 

Taking the solar system analogy to its extreme would suggest 
that the gap in the dust distribution seen around rj Tel at 5-20AU 
is caused by the existence of giant planets distributed in this 
region. Like the Solar System, the planetesimal belt at 3.9AU 
could represent a region where terrestrial planet formation was 
disrupted by the influence of a gas giant planet at 5-6AU, and 
the inner edge of the belt at 24AU could demark the outer edge 
of the planetary system. 

The inner edge of the the hot dust could also be evidence 
of terrestrial plan ets at <3AU, as the hot dust at lAU around 
HD 69830 ( Beichman et al.ll2005l:lLisse et al.l2007l) lies just out- 
side 3 Neptune m ass planets, the outermost of which is at 0.63 
AU (iLovis et al.ll2006i) (although the location of the disk has not 
yet been confirmed through resolved imaging). Also, although 
no planets have been detected around zeta Lep, the restricted 
spatial extent of its asteroid belt (Moerchen et al. 2007b) is sug- 
gested to be possible evidence of perturbations of planetary bod- 
ies. 

Implicit in the self-stirring interpretation is the formation of 
planets at least the size of Pluto in the inner regions. Such ob- 
jects would presumably have grown through mutual collisions 
to larger sizes, but it is not required that they also evolved into 
gas giant planets, which would depend on whether they rea ched 
~ lOMe before the gas disk dissipated (e.g. lLissaue 

A further possibility which is not discussed here is that this 
disk is a remnant of the protoplanetary disk, and that the dust is 
shepherded i nto a ring through the action of gas drag forces or 
instabilities (Takeu chi & Artvmowicj|2001l: iKlahr & Linll2005l; 
Besla & Wu 2007). In this case it would not be necessary for any 
planets to be present. A massive gas disk around a star of this 
age would be highly unus ual, as most gas disks a re believed to 
dissipate in < lOMyr (e.g. lHoUenbach et al.ll2005[). althou j gh gas 
has been detected around 4 A stars (J3 Pic'Brandeker et al. 2004^ 
HD 141596 Jonkheid et al. 2006; 49 Ceti Dent et al. 2005; and 
HD 32297 lRedfieldll2007i) . To assess this possibiHty knowledge 
of the gas distribution is vitally important, yet gas has yet to be 
detected around 77 Tel. 

5.5. Comparison with A star disks in /3 Pic moving group 

Of the 5 A stars in the f5 Pic moving group, 3 exhibit excess mid- 
IR emission from a debris disk (J3 Pic, 77 Tel and HD 172555). 
The derived (or directly measured) structures of these 3 disks 
display significant diversity. /3 Pic, an A6V star, has a well known 
highly extended disk with dust imaged from 10s of AU to an 
outer extent of ~ 1500AU ( T elesco et al.l2005l:lL^wood & KalasI 
2001). Since dust at the outer edge is th ought to be created i n 
planetesimal belts much closer to the star jAugereau et alj2001h . 
it is probably most instructive to consider this disk as an ex- 
tended planetesimal belt out to lOOAU with a bro ad peak in sur- 
face density at around 70AU dTelesco et al.ll2005l) . HD 172555 
on the other hand has all of its emission much closer to the star. 
Its emission spectrum indicates the presence of silicate grains at 
520K and a blackb ody com ponent at 200K im plying a radial 
offset of 0.8-6AU dChen et al. 2006; iRebuU e t al. 2008). The 
~24AU outer belt of the 77 Tel debris lies intermediate between 
these two extremes. These disks thus represent the diversity that 
is see n amongst debris d isks at the age when emission is seen to 
peak jCurrie et al.ll2007h . 

Within the context of the self-stirring model this diversity 
would be expected to arise primarily from different initial disk 
masses. While 77 Tel is the product of a 0.7x MMSN disk, the 
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protoplanetary disk of HD 172555 (and of the two A stars which 
are not detected) is presumably of lower density, and that of /3 
Pic is of higher density to allow Plutos to form out to 70AU by 
12Myr It is also worth noting that the beta Pic disk shows an 
asymmetry that may be the result of a catastrophic collision at 
~50AU, such as might be expected from the self-stirring models 
at the epoch of Pluto-formation. S uch a diversi ty of protoplane- 
tary disk masses is expected (e.g.. [Andrews ^Williams 2005]), 
but the small number of stars in this sample does not permit any 
definite conclusions to be drawn. In the young Solar System in- 
terpretation, the diversity in these disks could instead be due to 
the diversity in planetary system architecture across the three 
stars, which in add ition to a dependence on initial disk mass 
jWvatt et al.l2007a|) could also have a s ignificant stochastic com- 
ponent (e.g.^ Ouintana & Lissauej2006 , find small changes in the 
location of planetary embryos has a great impact on the final sys- 
tem architecture). Furthermore the origin of the hot emission, 
and its diversity, remains a puzzle, and it is not clear whether 
the presence, or lack thereof, of asteroid belts and ongoing ter- 
restrial planet formation should b e a monotonic function of the 
initial disk mass (see e.g. lSu et"al . 2008, who find from a sample 
of middle-aged disks around A0-A3V stars that factors beyond 
variation in initial disk mass are needed to explain the variation 
in observed disk structure). 

6. Conclusions 

We have presented Si-5 and Qa band Gemini T-ReCS imaging 
of the 12Myr old A star t] Tel. The emission at 18.3yL(m is shown 
to be significantly extended compared with a point source. This 
image represents the first resolution of dust emission around this 
star. These observations confirm the interpretation of the SED as 
a debris disk comprised of two populations: one component that 
is resolved at 24AU, the other that is unresolved with a temper- 
ature consistent with dust at 3.9AU. 

Two interpretations are proposed for the origin of the archi- 
tecture of these debris belts, (i) The system is a close analogue 
to the Solar System at a similar age, both in terms of the radial 
location of the debris and the level of emission; i.e., these belts 
represent analogous asteroid and Kuiper belts, ii) Alternatively 
this system is still undergoing planet formation, and the 24AU 
ring is the location where Pluto-sized objects recently formed; 
the 3.9AU ring could be the product of recent collisions between 
Mars-sized bodies in the final accumulation phases of terres- 
trial planets. This implies that this system evolved from a 0.7x 
MMSN protoplanetary disk. 

These possibilities can be further tested by probing the radial 
structure of the disk through multiple wavelength imaging, and 
by searching for planets both through direct methods, or through 
looking for evidence for them in the disk structure. Thus rj Tel is 
a valuable source for testing planet formation theories, not least 
since it may be an example of a system where planet formation 
is ongoing. 
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